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ABSTRACT 

A  probabilistic  event-step  computer  simulation  of  a 
repairable  item  inventory  system  with  exponential  inner- 
arrival  times  between  demands  is  constructed.   The  model 
allows  for  a  wide  variety  of  repair  times  and  lead  times 
to  be  considered.   The  major  parameters  are  investigated 
for  their  sensitivity  to  system  changes  in  order  to  help 
the  user  evaluate  given  or  proposed  inventory  policies 
and  parameters.   A  complete  description  of  the  model  is 
presented  and  various  measures  of  supply  performance  such 
as  the  expected  number  of  unit  years  of  backorders  per 
unit  period  and  the  fill  rate  are  provided.   All  flow 
charts  and  the  GPSS/360  program  listing  are  included. 
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I.   INTRODUCTION 

Traditional  inventory  analysis  is  ultimately  concerned 
with  two  decisions:  (1)  when  to  order  and  (2)  how  much  to 
order.  Most  of  the  work  done  in  the  past  addresses  itself 
to  these  questions  for  non-repairable  items  (consumables) . 
A  thorough  treatment  of  such  inventory  models  for  both  the 
backorders  and  lost  sales  case  for  numerous  resupply  and 
demand  distributions  now  exists. 

However,  at  least  in  military  applications,  repairable 
items  have  now  been  found  to  be  of  great  importance  both 
materially  and  financially.   Repairables  often  account  for 
a  significant  portion  of  the  items  controlled  by  an  activity 
and  even  more  significantly  when  their  portion  of  the  total 
dollar  investment  is  considered.   For  example,  focus  on  one 
of  the  inventory  control  points  of  the  U.S.  Navy  dealing 
with  repairables  such  as  the  Aviation  Supply  Office  (ASO) 
in  Philadelphia,  Pennsylvania.   At  fiscal  year's  end  1971 
repairables  represented  18.7%  of  the  different  types  of  line 
items  controlled  by  this  activity  for  a  total  of  66,000 
separate  items.   These  same  items  reflect  an  investment  of 
nearly  3.1  billion  dollars,  almost  60%  of  ASO's  total  dollar 
investment.   Further,  ASO  spends  about  16  million  dollars 
per  month  for  repairs  on  these  items  and  another  28  million 
dollars  monthly  for  procurement  of  such  items.   Among  ASO 
managed  repairable  items,  365,000  units  were  reworked  in 
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fiscal  year  1971.   In  view  of  the  above  figures  it  is  easy 
to  see  why  a  need  exists  to  help  the  manager  improve  his 
control  over  repairable  items.   It  appears  that  an  effective 
control  over  repairs  would  free  for  other  needs  many  ASO 
dollars  now  spent  on  procurement. 

Until  fairly  recently  the  control  of  repairable  items  had 
not  received  much  attention.   Even  in  the  case  of  consumables 
there  are  major  difficulties  in  applying  inventory  theory  to 
the  problem  because  of  the  intricate  relationships  of  the 
involved  item  parameters,  such  as  demand  and  lead  times. 
Also  the  number  of  different  types  of  line  items  complicates 
the  task.   The  introduction  of  repairable  items  makes  the 
problem  even  more  difficult  since  repairable  items  often 
require  special  storage,  handling,  packaging  and  shipping 
procedures  due  to  their  complexity.   In  addition,  the 
introduction  of  the  extra  source  of  supply,  the  repair 
department,  causes  further  complications.   On  the  other  hand, 
repairable  items  tend  to  minimize  the  number  of  different 
type  of  line  items  needed  to  be  carried  by  an  activity;  the 
parameters  are  usually  more  easily  determined  than  those  of 
consumables;  and  the  procurement  lead  times  are  more  reliable 
In  addition,  upper  bounds  for  the  total  demand  can  often  be 
determined,  for  if  a  given  repairable  item  is  a  subcomponent 
of  a  larger  unit,  the  number  of  larger  units  in  use  is 
generally  easy  to  determine. 

Few  people  have  addressed  the  problem  of  repairables. 
However,  Schrady  [Ref.  8]  and  others  at  the  U.S.  Naval 
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Postgraduate  School,  as  well  as  Feeney  and  Sherbrooke  [Ref.  6], 
have  looked  at  the  problem  for  the  deterministic  case. 
Feeney' s  and  Sherbrooke' s  model  included  the  capability  of 
handling  compound  Poisson  demands  but  without  degradation  of 
supply.   Wilson  and  Richards  [Ref.  10]  extended  this  model 
to  allow  for  this  possibility.   Due  to  problems  of  convexity 
and  others,  no  one  has  come  up  with  a  quick  and  easy  way  to 
evaluate  and  compare  alternative  repairable  item  inventory 
systems  on  a  larger  scale.   Richards  [Ref.  7]  addresses  this 
problem  and  this  paper  supplements  his  work  by  providing  a 
detailed  simulation  model  along  with  a  programming  guide  and 
flow  charts.   This  model  enables  the  user  faced  with  such 
problems  to  easily  come  up  with  numerical  values  for  various 
measures  of  performance  such  as  the  average  number  of  back- 
orders  per  unit  period  or  the  average  number  of  fills  per 
unit  period.   These  measures  of  performance  can  then  be  used 
by  the  manager  to  evaluate  and  compare  alternative  systems. 
Further,  the  model  can  also  be  used  as  a  teaching  aide  to 
help  students  and/or  supply  personnel  get  a  clearer  under- 
standing of  the  relationships  between  the  various  parameters 
involved  in  various  inventory  systems.   A  more  detailed 
description  of  the  model  and  how  to  use  it  are  given  in 
sections  II,  V,  and  VI  of  this  paper. 
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II.   THE  MODEL 

The  model  presented  in  this  paper  is  a  general  event-step 
computer  simulation  of  the  interactions  between  the  major 
parameters  for  Q  R,  continuous  review  inventory  systems. 
Here,  Q  is  the  order  quantity  and  R  is  the  reorder  point  for 
a  given  system.   The  purpose  of  this  general  model  is  to  test 
and  investigate  the  sensitivity  of  the  parameters  involved 
to  help  the  decision  maker  evaluate,  modify,  and/or  set  up 
appropriate  inventory  systems.   It  is  also  used  to  demonstrate 
the  analytical  results  obtained  by  Richards  [Ref.  7].   The 
inventory  systems  simulated  are  single  item,  single  echelon, 
and  repairable  or  non-repairable  systems  where  failures  are 
assumed  to  be  Poisson  distributed.   The  model  is  a  probabil- 
istic simulation  in  that  the  success  or  failure  of  any 
probabilistic  event  is  determined  in  the  model  by  comparing 
the  numerical  value  assigned  to  the  probability  of  success 
or  failure  to  a  program  generated  random  number.   In  this 
model  the  numbers  so  used  are  obtained  from  many  different 
distributions,  using  the  one  which  best  fits  the  particular 
parameter  being  considered  in  each  case.   Further,  it  is  an 
event-step  model  in  that  the  master  program  is  controlled  by 
a  subprogram  called  an  event  list  which  is  repeatedly  updated 
as  the  master  program  switches  from  one  event  subroutine  to 
another.   This  eliminates  the  necessity  to  look  at  the  status 
of  the  entire  system  after  each  time  unit  (one  hour)  to 
decide  whether  or  not  an  event  has  occurred. 
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Figure  1  shows  the  general  type  of  inventory  systems 
addressed  in  this  paper  .   In  order  to  help  the  reader  under- 
stand the  relationships  pictured  there,  the  following 
definitions  for  the  on-hand  inventory,  repair  inventory, 
and  others  are  given.   Note  that  because  of  the  addition  of 
the  repair  facility,  the  commonly  used  definitions  of  the 
various  states  for  the  consumable  item  inventory  system 
must  be  modified. 

Definition  2.1:   The  on-hand  inventory  (D-B)  includes 
all  of  the  items  ready  for  issue  (D)  less  any  backorders  (B) . 
Note  that  the  on-hand  inventory  will  be  positive  if  there  are 
no  items  on  backorder. 

Definition  2.2:   The  repair  inventory  (C)  includes  all 
items  at  the  repair  facility,  including  both  those  undergoing 
repair  and  those  waiting  to  begin  repair. 

Definition  2.3:   The  net  inventory  (D-B+C)  includes  the 
on-hand  inventory  and  the  repair  inventory. 

Definition  2.4:   The  inventory  position  (D-B+C+A)  is 
the  sum  of  the  net  inventory  plus  all  items  which  have  been 
placed  on  order  but  have  not  yet  been  received. 

Consider  the  inventory  system  as  one  which  has  control 
of  a  large  number  of  items  of  a  single  type  where  the 
number  of  demands  (failures)  for  any  time  period  is  a  value 
of  a  random  variable  having  a  Poisson  distribution  with 
mean  rate  %.      The  items  which  fail  are  assumed  to  be 
repairable  with  probability  P  and  with  probability  (1-P) 
the  items  are  scrapped.   These  failures  will  be  referred  to 
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Figure  I.    Repairable  item  inventory  system. 
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as  type  one  and  type  two,  respectively.   The  manager 
controlling  the  system  attempts  to  minimize  the  time 
required  to  fill  a  demand  by  maintaining  a  stock  of  items 
ready  for  issue.   Systematically,  the  manager  supplements 
his  stock  by  procurement  of  new  items  and/or  receipt  of 
repaired  items.   Hence  the  manager  is  faced  with  the  problems 
of  when  and  how  much  to  reorder.   Standard  approaches  try 
to  answer  these  questions  by  concentrating  their  efforts  on 
minimizing  the  total  expected  cost.   However,  a  difficulty 
in  this  approach  is  the  determination  of  proper  weights  to 
assign  to  the  various  cost  factors  involved. 

A  second  approach  concentrates  on  a  class  of  simple 
procurement  policies  based  on  one  or  two  parameters.   The 
goal  of  this  approach  is  to  describe  the  long-term  effect 
of  an  arbitrary  member  of  the  class  of  simple  procurement 
policies  being  considered  and  then  determine  the  choices 
for  the  corresponding  parameters  which  are  best  with  respect 
to  a  given  measure  of  performance.   This  computer  simulation 
allows  the  user  to  do  this  quickly  and  easily. 

After  a  class  of  simple  policies  has  been  specified, 
the  study  of  the  inventory  model  becomes  a  study  of  the 
associated  stochastic  processes  and  the  long  run  properties 
of  the  processes.    Emphasis  is  given  to  determination  of 
stationary  distributions  which  describe  the  probabilities 
of  being  in  the  various  states  after  the  system  has  been  in 
operation  for  a  long  period  of  time.   This  computer  program 
greatly  simplifies  the  procedure  whereby  actual  numbers  are 
produced  for  these  various  entities. 
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The  program  is  capable  of  representing  the  repair 
department  in  one  of  two  ways.   First,  as  an  infinite- 
server  queue  with  general  independent  repair  times  which 
can  be  constant  or  realizations  of  a  uniform,  normal,  or 
empirical  distribution.   Under  this  policy  the  repair 
department  acts  as  a  M/G/«©  queueing  system  which  for  the 
random  repair  times  case  will  allow  crossing  of  repair  times 
The  crossing  of  orders  is  also  allowed  in  this  case. 

The  second  method  considers  the  repair  department  as  a 
single  server  queueing  system.   That  is,  only  one  batch  of 
items  up  to  a  maximum  quanitity  (K)  can  be  repaired  at  one 
time.   Under  this  policy  the  user  also  has  the  option  to 
enforce  a  constraint  on  the  minimum  batch  size  (L)  as-  well. 
This  latter  policy  is  provided  because  frequently  in  the 
case  of  repair,  large  set-up  costs,  tool-up  costs,  or 
manpower  constraints  make  the  second  repair  policy  more 
economical. 
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III.   EXPERIMENTS 

A.   OBJECTIVES  AND  APPROACH 

1.  Objectives 

The  objectives  of  this  paper  are  threefold.   (1)  To 
show  the  need  for  and  justification  of  such  a  simulation. 
(2)  To  demonstrate  the  analytical  results  obtained  by 
Richards  [Ref .  7] .   (3)  To  provide  a  means  through  which  a 
decision  maker  can  evaluate  given  or  hypothetical  inventory 
systems. 

The  first  objective  has  been  touched  on  already  and 
will  be  complemented  by  the  rest  of  this  paper.   Consistent 
with  the  second  and  third  objectives,  a  number  of  experiments 
were  performed  which  involved  generating  operating  statistics 
in  order  to  quantitatively  compare  the  simulated  results  to 
the  theoretical  and  display  the  sensitivity  of  various 
parameters  for  particular  systems. 

2.  Approach 

Each  experiment  consisted  of  a  series  of  simulation 
runs.   Each  run  simulated  up  to  14  unit  periods,  where  a 
unit  period  represented  ten  years  of  actual  operation  of  a 
given  inventory  system.   For  each  run,  operating  statistics 
were  collected  for  all  measures  of  performance  listed  in 
Table  I.   For  illustration  purposes  two  measures  of  per- 
formance were  chosen  to  judge  the  performance  of  the  hypo- 
thetical systems  examined  in  this  paper.   Those  chosen  were 


20 


TABLE  I.   SUPPLY  PERFORMANCE  MEASURES 


MOP   1   (F)      =  XH100:   The  average  number  of  fills  in  a 
unit  period. 

MOP   2   (FILRT)   =  XH99  :   The  fill  rate  multiplied  by  100. 
This  is  the  probability  of  being  able  to  satisfy  a 
demand  without  a  time  delay. 

MOP   3   (B)      =  XH98  :   The  average  number  of  backorders 
incurred  in  a  unit  period. 

MOP   4   (B^)     =  XH97  :   The  expected  number  of  unit  years 
of  backorders  in  a  unit  period  multiplied  by  100. 

MOP   5   (BARB)    =  XH96  :   The  average  number  of  items  on 
backorder  at  an  arbitrary  point  in  time  multiplied 
by  100. 

MOP   6   (RDYRT)   =  XH95  :   The  probability  of  being  able  to 
fill  a  demand  at  an  arbitrary  point  in  time  multi- 
plied by  100. 

MOP   7   (ENOH)    =  XH94  :   The  average  number  of  items  on 

hand  at  an  arbitrary  point  in  time  multiplied  by  100, 

MOP  8  (ORD)  =  XH93  :  The  average  number  of  orders 
placed  in  a  unit  period. 

MOP  9  (RI)  =  XH92  :  The  average  number  of  repair 
inductions  (batches  of  items  repaired)  per  unit 
period.  In  the  infinite  server  queueing  system 
RI  will  be  C. 

MOP  10   (C)      =  XH91  :   The  average  number  of  items  in 

the  repair  department  at  an  arbitrary  point  in  time 
multiplied  by  100. 

MOP  11   (DIFF  1)  =  XH90  :   The  difference  between  ENOH  and 

the  theoretical  value  (V4)  of  the  expected  number  of 
items  on  hand  at  an  arbitrary  point  in  time  multi- 
plied by  100. 

MOP  12   (DIFF  2)  =  XH89  :   The  difference  between  BT  and 
BARB  multiplied  by  100.   This  is  calculated  to 
illustrate  that  the  total  time-weighted  backorders 
incurred  in  a  unit  period  is  equivalent  to  the 
expected  number  of  backorders  at  an  arbitrary  point 
in  time. 
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NOTE:   Measure  of  performance  four  (BT)  would  normally  be 
normalized  (TTT)  by  dividing  by  ten,  i.e., 
XH97  =  1102=  ($£)  =1.102. 

the  expected  number  of  unit  years  of  backorders  and  the 
average  number  of  items  on  hand  at  an  arbitrary  point  in 
time.   For  all  runs  the  system  was  preloaded  such  that  the 
inventory  position  was  equal  to  the  sum  of  the  order  quantity 
and  reorder  point.   It  was  felt  that  this  would  result  in 
negligible  transient  times,  hence  the  entire  run's  statistics 
were  usable  in  all  cases.   The  input  values  of  all  parameters, 
for  what  was  called  the  parent  system,  are  summarized  in 
Table  II.   Deviations  from  these  values  of  the  parent 
system  for  any  one  experiment  are  pointed  out  in  the 
description  of  each  experiment  as  well  as  in  Table  III. 

TABLE  II.   PARENT  SYSTEM  PARAMETER  VALUES 
PARAMETER  VALUE 

Order  quantity  (Q)  =10 

Reorder  point  (R)  =0,1,2,3... 

Minimum  batch  size  (L)  =0 

Maximum  batch  size  (K)  =0 

Probability  of  repair  (P)  =.9 

Mean  demand  during  repair  time  (X)        =1.4 
Type  of  distribution  for  repair  time      -constant 
Mean  repair  time  (X  )  =336  hours 

Lower  bound  for  uniform  repair  time  (LBR)  =0 
Upper  bound  for  uniform  repair  time  (UBR)  =0 
Variance  for  repair  time  distribution  (CR)=0 
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Type  of  distribution  for  lead  time  -constant 

Mean  lead  time  (X  )  =504 

Li 

Lower  bound  for  uniform  lead  time  (LBL)  =0 

Upper  bound  for  uniform  lead  time  (UBR)  =0 

Variance  for  lead  time  distribution  (6^)  =0 

Type  of  output  desired  -summary 

TABLE  III.   CHANGES  TO  PARENT  SYSTEM  PARAMETER  VALUES 


A.   CHANGES  FOR  EXPERIMENT  ONE 


B.   CHANGES  FOR  EXPERIMENT  TWO 


(X)=1.4,  2.75,  4,  5.5 


Type  of  distribution  for  repair  time--constant , 

uniform,  normal, 
and  empirical 

Type  of  distribution  for  lead  time   --constant, 

uniform,  normal, 
and  empirical 

LBR=168  hours 

UBR=504  hours 

6r  =  75  hours 

C.  CHANGES  FOR  EXPERIMENT  THREE 
(Q)=2  (K)=l,  2,  3 

D.  CHANGES  FOR  EXPERIMENT  FOUR 
(Q)=2  (L)=l,  2,  3,  4 

E.  CHANGES  FOR  EXPERIMENT  FIVE 
(QJ=2  (K)=4 
(T)=0.40,  0.65,  0.75,  0.80,  0.90,  1.15,  1.40 


LBL=240 

UBL=768 

ETl  =116 

,    5,    6 

(L)=l 

(10=4 

(L)=l 
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B.   DESCRIPTION  OF  EXPERIMENTS 

1.   First  Experiment:   The  Whole  Spectrum 

The  first  experiment  considered  the  parent  system 
but  with  values  of  mean  demand  during  repair  time  of  1.4, 
2.75,  4,  and  5.5.   Also,  order  quantities  of  one  through 
ten  and  variable  reorder  points  were  allowed.   That  is,  for 
each  combination  of  demand  and  order  quantity  above,  the 
reorder  point  (R)  was  allowed  to  vary  from  an  initial  value 
of  zero  until  the  expected  number  of  unit  years  of  backorders 
reached  an  acceptable  stationary  level.   This  was  done  in 
order  to  see  how  the  measures  of  performance,  particularly 
the  two  mentioned  above,  varied  with  respect  to  changes  in 
the  order  quantity  and  reorder  point  for  the  system  described 
above. 

The  amount  of  computer  time  required  to  run  one  period 
ranged  from  a  low  value  of  5.89  seconds  for  (X)=1.4,  with 
Q=10  and  R=7;  to  a  high  value  of  28.57  seconds  for  (X)=5.5, 
with  Q=l,  and  R=0.   The  total  computer  time  required  for 
this  experiment  would  obviously  depend  on  the  particular 
inventory  system  being  investigated  by  the  user.   It  is 
assumed  that  the  user  would  know,  determine,  or  estimate 
his  system  demand  such  that  he  could  restrict  his  investiga- 
tion to  one  or  two  values.   Therefore,  the  number  of  runs 
and  total  computer  time  required  for  his  experiments  would 
be  significantly  less  than  that  required  for  these  presented 
in  this  paper.   Some  typical  results  of  the  first  experiment 
are  shown  in  Figures  2  and  3. 
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2.   Second  Experiment:   Distribution  Test 

For  this  experiment  a  representative  system  from 
experiment  one  was  chosen  for  further  study.   The  system 
chosen  was  thought  to  be  a  typical,  low  demand,  high  cost, 
repairable  item  inventory  system  with  a  reorder  quantity  of 
ten  and  a  mean  demand  during  repair  time  of  5.5.   The  par- 
ticular reorder  quantity  used  is  a  somewhat  subjective 
choice  which  the  manager  would  make.   Hopefully  the  manager 
could  look  at  the  results  of  experiment  one  and  pick  what 
would  be  the  best  quantity.   A  procedure  to  do  this  is 
explained  in  the  analysis  of  experiment  one.   However,  what 
will  often  happen  is  that  one  of  several  possible  factors 
will  force  the  manager  to  live  with  some  predetermined  order 
level.   Some  examples  of  these  may  be:   there  is  a  minimum 
number  of  items  that  the  manager  can  let  a  contract  for  with 
a  particular  manufacturer,  or  an  engineering  requirement 
forces  items  to  be  procured  in  matched  pairs  or  sets,  or 
possibly  a  fixed  unit  of  issue  by  the  manufacturer.   Regard- 
less of  how  it  is  decided  upon,  once  the  reorder  quantity  is 
fixed  the  sensitivity  of  the  system  can  be  investigated  by 
changing  various  input  factors  such  as  the  repair  time 
distribution  or  lead  time  distribution,  as  was  done  in  this 
experiment.   Constant,  uniform,  normal,  and  empirical  repair 
and  lead  time  distributions  were  examined.   For  each  combina- 
tion of  distributions  the  reorder  point  (R)  was  varied  until 
the  total  unit  years  of  backorders  was  stationary.   The 
purpose  of  this  experiment  was  to  show  the  sensitivity  or 
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insensitivity  of  the  modified  parent  system  to  various 
types  of  distributions  for  the  repair  and  lead  times  each 
having  the  same  mean  and  variance.   Some  typical  results  of 
the  experiment  are  given  in  Figure  4. 

3.   Third  experiment:   Maximum  Batch  Size  Test 

For  this  experiment  the  parent  system  was  modified 
in  order  to  treat  the  repair  department  as  a  single  server 
queueing  system  where  the  server  is  capable  of  handling  as 
many  as  K  items  at  one  time.   Using  the  procedure  explained 
in  the  analysis  of  experiment  one,  the  order  quantity  was 
fixed  at  a  value  of  two  for  this  experiment.   Then  the  sen- 
sitivity of  the  reorder  point  (R)  was  investigated  by  varying 
R  for  each  feasible  value  of  K,  the  maximum  batch  size,  until 
the  expected  number  of  unit  years  of  backorders  reached  an 
acceptable  stationary  level.   A  simple  way  to  determine  the 
maximum  value  for  K  is  to  look  at  the  results  obtained  in 
experiment  one  for  the  comparable  infinite  server  queueing 
system.   The  run  in  experiment  one  where  all  the  input 
parameter  values  were  identical  to  this  one  with  the  excep- 
tion of  the  parameters  associated  with  the  repair  discipline 
would  be  the  comparable  infinite  server  queueing  system.   One 
of  the  items  included  in  the  summary  output  for  Storage 
Statistics,  FASTR,  summarizes  the  statistics  for  the  case 
in  which  the  repair  department  is  treated  as  an  infinite 
server  queue.   Two  of  the  particular  statistics  included 
there  are  the  average  number  of  items  in  the  repair  depart- 
ment and  the  maximum  number  of  items  that  were  ever  in  the 
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repair  department  at  any  time.   These  particular  values  were 
used  in  this  experiment  as  a  guide  for  determining  the 
maximum  value  for  K  to  be  considered.   The  results  of  this 
experiment  are  shown  in  Figure  5. 

4.  Fourth  experiment:   Minimum  Batch  Size  Test 

In  this  experiment  the  parent  system  was  again 
modified  to  treat  the  repair  department  as  a  single  server 
queueing  system.   Again,  as  in  the  last  experiment,  the 
sensitivity  of  the  reorder  point  to  the  type  of  repair 
policy  used  was  investigated.   This  time,  however,  K 
rather  than  L  was  held  constant  while  R  was  varied  for  each 
value  of  L,  the  minimum  batch  size,  until  the  expected 
number  of  unit  years  of  backorders  reached  an  acceptable 
stationary  level.   The  values  of  L  were  varied  from  one  to 
K.   The  value  of  K  used  was  four,  which  was  the  maximum 
feasible  quantity  for  K  as  determined  by  experiment  three. 
Since  the  repair  department  was  treated  as  a  single  server 
queueing  system,  the  repair  of  one  batch  of  items  must  be 
completed  before  repairs  can  begin  on  the  next  batch  of 
items.   The  results  of  this  experiment  are  presented  in 
Figure  6. 

5.  Fifth  Experiment:   Mean  Repair  Time  Test 

The  modified  parent  system  of  experiment  four  was 
used  for  this  experiment  but  with  L=l  and  with  various  mean 
demands  during  repair  time  (X) .   The  mean  demand  during 
repair  time  was  allowed  to  vary  from  1.40  to  0.40  by 
changing  the  mean  time  to  repair  from  336  hours  to  96  hours 
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while  keeping  the  mean  time  between  failures  constant  at 
242  hours.   The  sensitivity  of  the  reorder  point  was  then 
investigated  by  again  varying  R  for  each  value  of  X  until 
the  expected  number  of  unit  years  of  backorders  per  unit 
period  reached  a  stationary  level.   A  physical  interpretation 
here  could  be  a  change  of  equipment  (perhaps  manual  repair  to 
automation)  or  a  change  in  the  budget  for  the  repair  depart- 
ment.  An  investigation  of  the  results  of  this  experiment 
suggested  a  comparison  of  the  trade  off  between  decreasing 

the  mean  repair  time  (X  )  and  increasing  the  maximum  batch 

R 

size  in  the  repair  department.   The  results  of  experiment 
five  are  presented  in  Figure  7. 
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IV.   ANALYSIS  AND  CONCLUSIONS 

A.   ANALYSIS  OF  THE  FIRST  EXPERIMENT 

In  all  experiments  the  reorder  quantity  and  reorder 
level  were  controlled.   In  addition,  for  each  experiment 
one  variable  such  as  the  maximum  or  minimum  repair  batch 
size  was  controlled  while  the  remaining  variables  were 
determined  from  the  simulation.   In  this  experiment,  for 
each  reorder  quantity  and  mean  demand  during  repair  time 
considered,  the  reorder  level  which  reduced  the  expected 
number  of  unit  years  of  backorders  per  unit  period  to  an 
acceptable  level  was  determined.   To  determine  this  reorder 
point  a  one-way  analysis  of  variance  was  used. 

Such  a  test  requires  additivity,  linearity,  normality, 
independence,  and  homogeneous  variances.   The  first  two 
assumptions  were  felt  to  be  satisfied  since  the  expected 
number  of  unit  years  of  backorders  per  unit  period  appears 
to  be  a  linear  sum  of  the  true  mean  effect,  the  true  effect 
of  the  treatment  (the  value  of  the  reorder  point  chosen) , 
and  the  true  effect  of  the  errors  observed.   The  primary 
error  in  the  observations  is  due  to  the  stream  of  random 
numbers.   It  seems  to  be  a  good  assumption  that  these  errors 
are  normally  distributed  with  mean  zero  and  variances  which 
are  equal  for  each  treatment. 

This  particular  test  was  used  because  it  is  quite  robust, 
that  is,  the  test  is  reliable  under  rather  strong  modifications 
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to  the  above  assumptions.   It  can  also  be  shown  that  trans- 
formations of  the  original  data  are  justified.   This  was 
particularly  important  in  this  case  because  it  was  necessary 
to  make  a  logarithmic  transformation  to  the  observations 
in  order  to  obtain  a  linear  model.   Figure  8  illustrates 
this  type  of  transformation. 

The  null  hypothesis  used  in  this  analysis  of  variance 
was  that  the  expected  number  of  unit  years  of  backorders 
per  unit  period  for  each  of  the  reorder  points  was  the  same 
for  a  fixed  reorder  quantity  and  mean  demand  during  repair 
time.   For  the  parameters  considered  in  this  paper  the 
analysis  of  variance  was  reported  for  that  case  which 
appeared  to  have  the  highest  probability  of  not  rejecting 
the  null  hypothesis.   Nine  observations  with  different 
streams  of  random  numbers  for  each  run  were  made  and  the 
null  hypothesis  was  tested  at  the  .01  significance  level. 
The  results  of  the  analysis  of  variance  yielded  a  test 
statistic  of  194,  well  in  excess  of  the  F-value  (2.48) 
required  to  reject  the  null  hypothesis.   (For  calculations 
see  page  79.) 

Since  the  primary  goal  of  the  first  experiment  was  to 
determine  the  "best"  order  quantity,  reorder  point  combina- 
tion, Scheffe's  [Ref.  5]  multiple  comparison  test  was  then 
used  to  find  this  combination.   This  procedure  tests  for 
equality  of  means  between  contrasts  in  a  systematic  fashion 
until  a  significant  difference  between  means  is  found.   The 
analysis  of  variance  indicated  that  a  significant  difference 
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was  present.   Scheffe's  test  was  then  used  to  determine 
the  minimum  value  of  the  difference  in  means  which  statistic- 
ally differentiated  among  the  various  values  of  the  expected 
number  of  unit  years  of  backorder  per  unit  period. 

In  order  to  apply  Scheffe's  test  a  (Q,R)  pair  yielding  a 
value  of  the  measure  of  performance  which  was  considered 
acceptable  was  chosen.   The  first  contrast  looked  at  the 
differences  in  the  measure  of  performance  for  this  pair 
(Q,R^)  and  the  pair  (Q,RM~1) .   If  this  contrast  failed  to 
reject  the  null  hypothesis,  that  the  expected  number  of 
unit  years  of  backorders  per  unit  period  were  equal  at  the 
.01  significance  level,  then  the  next  contrast  would  examine 
the  pairs  (Q,RM)  and  (Q,Rw-2).   This  procedure  continued 
until  the  null  hypothesis  was  rejected.   If  this  happened 
when  the  pairs  (Q,R^)  and  (Q,R  -j)  were  compared,  the 
minimum  value  of  the  differences  which  was  statistically 
significant  was  known  to  be  less  than  or  equal  to  the 
difference  between  the  values  of  the  measure  of  performance 
for  these  two  pairs.   For  future  reference  denote  Rw-j  by  R 
The  minimum  difference  required  for  significance  can  be 
solved  for  exactly  by  applying  this  procedure  in  reverse, 
that  is,  by  setting  the  value  of  an  arbitrary  contrast 
equal  to  the  test  statistic  required  to  reject  the  null 
hypothesis  and  then  transforming  this  contrast  into  the 
difference  required  in  the  measure  of  performance  for 
significance.   For  the  example  systems  considered  in  this 
paper  this  procedure  required  the  use  of  three  contrasts 
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before  significance  was  found  between  (Q,R)  pairs  (10,11) 
and  (10,8).   This  contrast  resulted  in  a  test  statistic  of 
6.99  which  was  greater  than  the  F-value  (2.48).   Hence,  the 
null  hypothesis  was  rejected.   (For  calculations  see  pages 
81-83).   Further,  forcing  the  test  statistic  to  be  equal  to 
2.48  and  solving  for  the  required  minimum  difference  between 
values  of  the  measure  of  performance  for  significance 
resulted  in  a  value  of  0.096  unit  years  of  backorders  per 
unit  period.   Statistically,  one  cannot  distinguish  a 
difference  between  0.096  and  0.00  for  this  measure  of 
performance.   For  this  particular  example,  which  is  summar- 
ized in  Tables  IV  and  V,  this  means  that  any  position  to 
the  right  of  the  dashed  line  in  Table  IV  could  be  an 
acceptable  operating  point  for  a  system  when  X,  the  mean 
demand  during  repair  time,  was  5.5  and  the  other  system 
characteristics  were  as  noted  in  Tables  II  and  III. 

The  procedure  then  used  to  determine  the  "best"  operating 
point  for  this  particular  example  was  to  list  the  pairs 
(Q,R  )  and  observe  the  measures  of  performance  listed  in 
Table  I.   Notice  from  Table  V  that  when  one  tries  to 
identify  the  pair  that  minimizes  both  the  average  number  of 
items  on  hand  at  an  arbitrary  point  in  time  and  the  average 
number  of  procurements  per  unit  period  the  (Q,R0)  pairs 
(2,10)  and  (3,10)  are  selected  as  the  "best"  operating 
points  for  this  example.   A  simple  comparison  of  holding 
cost  and  procurement  cost  could  enable  the  user  to  further 
distinguish  between  these  two  choices. 
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In  applying  these  results  to  other  systems  it  was  ob- 
served that  the  value  of  0.096  for  the  expected  number  of 
unit  years  of  backorders  per  unit  period  could  also  be  used 
to  determine  the  (Q,R„)  candidates  for  the  "best"  operating 
point  when  considering  all  values  of  mean  demand  during 
repair  time  less  than  or  equal  to  5 . 5 .   It  was  further 
observed  that  among  these  (Q,R  )  candidates,  most  measures 
of  performance  would  differ  by  negligible  amounts  with  the 
exception  of  the  average  number  of  items  on  hand  at  an 
arbitrary  point  in  time  and  the  average  number  of  procure- 
ments made  per  unit  period. 

B.   ANALYSIS  OF  EXPERIMENTS  TWO  THROUGH  FIVE 

1.  General  Comments 

The  analysis  of  variance  and  the  multiple  comparison 
test  could  have  been  used  to  obtain  candidates  for  optimal 
operating  points  with  other  relevant  parameters,  such  as 
the  minimum  or  maximum  repair  batch  sizes  or  the  distribution 
of  repair  time  being  used  as  the  basis.   Although  this 
analysis  was  not  repeated  for  all  experiments,  a  close 
observation  of  the  data  seemed  to  indicate  that  the  minimum 
significant  difference  value  would  be  nearly  the  same 
regardless  which  parameter  was  chosen  for  the  basis.   Use 
of  this  observation  was  made  in  the  analysis  of  experiments 
two  through  five. 

2 .  Repair  Time  Distribution  (Second  Experiment) 

The  results  summarized  in  Tables  VI  and  VII  showed 
that  the  type  of  repair  time  distribution  had  a  negligible 
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effect  on  the  expected  number  of  unit  years  of  backorders  per 
unit  period  provided  the  mean  repair  times  were  the  same. 
This  result  was  expected  since  Palm's  theorem  indicates  that 
the  distribution  of  the  number  of  items  in  repair  depends 
only  on  the  repair  time  distribution  through  its  mean  value. 
When  the  probability  of  repair  was  nearly  one  it  was  also 
observed  that  the  type  of  lead  time  distribution  had  a 
negligible  effect  on  system  performance.   A  probable 
explanation  for  this  result  is  that  the  average  number  of 
orders  placed  per  unit  period  was  insignificant  compared  to 
the  number  of  repairs  made.   This  can  be  seen  in  Table  VII. 
3.   Maximum  Batch  Size  Test  (Third  Experiment) 

In  this  experiment  the  effects  of  change  in  the 
maximum  repair  batch  size  were  investigated.   The  results 
summarized  in  Tables  VIII  and  IX  indicate  a  vast  improvement 
in  system  performance  for  increasing  the  maximum  batch  size 
from  one  to  two  but  negligible  improvement  thereafter.   The 
initial  improvement  can  be  explained  by  the  queueing  result 
which  says  that  in  order  to  ensure  stationarity  the  maximum 
repair  batch  size  must  be  large  enough  to  guarantee  that  the 
mean  time  between  failures  multiplied  by  K  must  be  greater 
than  the  mean  repair  time.   An  important  question  which  was 
considered  was  exactly  how  large  K  should  be.   In  order  to 
provide  an  answer  further  work  would  have  to  be  done  to 
determine  an  explicit  rule  for  finding  the  "best"  value  of  K. 
However,  it  is  felt  that  this  rule  would  depend  on  the 
utilization  factor  of  the  repair  department.   For  this 
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particular  example,  increasing  K  to  obtain  a  utilization 
factor  of  less  than  .417  resulted  in  negligible  improvement 
in  system  performance. 

4.  Minimum  Batch  Size  Test  (Fourth  Experiment) 

This  experiment  considered  the  effects  of  change  in 
the  minimum  batch  size  for  a  fixed  value  of  K.   The  results 
summarized  in  Tables  X  and  XI  indicated  that,  as  before, 
changes  in  the  minimum  batch  size  requirement  had  little 
effect  in  overall  system  performance  with  the  exception  of 
the  average  number  of  repair  inductions.   If  the  cost  of  a 
repair  induction  is  large  certain  economies  could  be  achieved 
by  requiring  a  minimum  batch  size.   This  could  be  achieved 
without  degrading  system  performance.   Further  investigation 
would  be  needed  to  generalize  these  results  in  order  to 
obtain  a  rule  for  determining  the  "best"  value  for  the 
minimum  batch  size. 

5.  Mean  Repair  Time  Test  (Fifth  Experiment) 

For  the  case  in  which  both  K  and  L  are  one,  the 
effects  of  change  in  the  mean  repair  time  were  examined. 
The  results  summarized  in  Tables  XII  and  XIII  indicated  that, 
as  expected,  the  magnitude  of  the  mean  repair  time  is  highly 
significant.   In  particular,  system  effectiveness  was  inversely 
proportional  to  the  mean  repair  time. 
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C.   CONCLUDING  COMMENTS  ABOUT  THE  ANALYSIS 

It  is  interesting  to  note  the  tradeoffs  between  various 
parameters.   For  instance,  it  seemed  that  one  could  deter- 
mine a  reduction  in  mean  repair  time  which  would  give  an 
equivalent  increase  in  system  performance  to  that  provided 
by  an  increase  in  the  maximum  repair  batch  size.   An  illus- 
tration of  such  a  tradeoff  for  the  example  considered  in 
this  paper  can  be  seen  by  looking  at  Table  XIII  for  the 
following  entries:   minimum  and  maximum  batch  size  of  one 
with  mean  repair  time  of  216  versus  minimum  batch  size  of 
one,  maximum  batch  size  of  four,  and  mean  repair  time  of 
336  hours. 

Decisions  such  as  determining  the  maximum  batch  size  or 
mean  repair  time  must  be  made  by  the  decision  maker  based  on 
the  feasibility  of  change  and  the  cost  of  such  a  change. 
This  simulation  will  provide  the  manager  with  each  of  the 
measures  of  performance  needed  to  make  such  decisions. 

In  addition  to  the  above  remarks,  another  item  should  be 
discussed.   Recall  that  one  of  the  objectives  of  this  study 
was  to  demonstrate  the  analytical  results  obtained  by 
Richards  [Ref .  7] .   Notice  that  in  each  of  the  tables  which 
summarized  the  measures  of  performance  for  the  various 
experiments,  two  of  the  measures  of  performance  listed  are 
"DIFF  1"  and  "DIFF  2."   "DIFF  1"  reflects  the  difference 
between  the  simulated  and  theoretical  values  for  the  average 
number  of  items  on  hand  at  an  arbitrary  point  in  time. 
"DIFF  2"  reflects  the  difference  between  the  expected  number 


50 


of  unit  years  of  backorders  per  unit  period  and  the  average 
number  of  items  on  backorder  at  an  arbitrary  point  in  time. 
In  observing  these  differences  in  Table  V  it  was  noted  that 
the  maximum  value  for  "DIFF  1"  was   0.39  and  for  MDIFF  2" 
it  was  .009.   This  represents  a  maximum  percentage  differ- 
ence of  5.1%  for  "DIFF  1."   It  is  felt  that  since  these 
differences  arc  so  small,  this  simulation  does  in  fact 
demonstrate  the  validity  of  the  analytical  results  obtained 
by  Richards. 

V.   PROGRAM  GUIDE 

This  section  is  written  so  that,  along  with  Table  I 
(Supply  Performance  Measures)  and  the  table  of  symbols  and 
abbreviations,  it  can  be  a  self  contained  unit  for  the  user, 
It  contains  all  the  information  needed  to  successfully  set 
up  and  run  the  model. 

A.   GENERAL  COMMENTS 

Most  people  will  agree  that  any  inventory  system  should 
be  designed  to  give  the  maximum  amount  of  service  for  the 
minimum  amount  of  cost  possible.   Hence  most  inventory 
studies  evaluate  the  various  policies  being  considered  by 
looking  at  the  costs  incurred  over  a  unit  period  of  time. 
Then  policy  A  is  chosen  over  policy  B  if  the  total  expected 
cost  per  unit  time  using  policy  A  is  less  than  when  using 
policy  B.   However,  quite  often  in  this  approach  one  cannot 
specify  all  the  cost  elements.   An  alternative  approach  is 
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to  determine  the  average  values  of  such  measures  of  perform- 
ance as  the  number  of  items  on  hand  at  an  arbitrary  point 
in  time,  the  number  of  backorders  per  unit  period,  and  others 
as  listed  in  Table  I.   The  user  can  then  choose  from  these 
measures  of  performance  those  he  wishes  to  use  to  evaluate 
the  policies  under  his  scrutiny.   This  method  provides  the 
user  with  an  easy  way  of  deciding  which  values  for  system 
parameters  such  as  the  reorder  quantity  and  reorder  point 
yield  satisfactory  performance. 

Closed  form  solutions  can  be  obtained  for  the  various 
measures  of  performance  but  often  the  solutions  are  not  very 
tractable  for  further  calculations.   This  is  where  the 
simulation  is  of  value.   It  systematically  and  quickly 
calculates  the  magnitude  of  all  the  measures  of  performance 
for  various  input  values  of  all  parameters  involved  and 
presents  the  user  with  a  fast  and  easy  means  of  investigating 
the  sensitivity  of  his  inventory  system  to  any  changes  in 
these  parameters.   To  try  and  do  this  procedure  by  hand  or 
apply  the  analytical  results  to  a  sensitivity  analysis  of 
this  size  would  involve  a  large  number  of  man-hours  of  work 
to  produce  effective  results. 

B.   EVENTS 

This  model  is  an  event-step  computer  simulation  written 
in  the  GPSS/360  language,  i.e.,  all  actions  that  are  to 
occur  in  the  simulation  are  dynamically  generated  by  the 
computer  program  as  a  result  of  previous  simulation  actions 
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and  are  listed  chronologically  in  order  of  descending 
priority  in  a  Current  Events  Chain,  Future  Events  Chain, 
or  Interrupt  Chain.   For  purposes  of  programming,  each  of 
the  major  actions  included  in  the  simulation  can  be  con- 
sidered as  an  entity  in  itself  and  assumes  the  form  of  a 
computer  program  subroutine.   Such  major  actions  have  been 
called  events  for  convenience.   The  eight  major  actions 
(events)  included  in  the  simulation  are: 

(1)  Item  Failure 

(2)  Item  Inspection 

(3)  Scrap 

(4)  Repair 

(5)  Issue 

(6)  Backorder 

(7)  Order 

(8)  Receipt  of  Material 

Each  of  the  computer  program  subroutines  representing 
these  events  uses  as  input  parameters  the  following  informa- 
tion: 

(1)  Time  the  event  is  to  occur 

(2)  Identification  of  event 

(3)  Number  of  items  involved 

(4)  Priority  level  (0,1,...) 

(5)  Current  values  of  system  characteristics  such  as 
inventory  position  (IP) ,  number  of  items  ready  for 
issue  (D) ,  and  others. 

The  dynamic  process  of  simulating  one  inventory  policy 

from  start  to  finish  forms  the  executive  routine  for  the 

computer  simulation.   This  executive  routine  also  includes 

three  program  subroutines  referred  to  as  Random  Observation, 

Calculation  of  Measures  of  Performance,  and  Incrementation 

of  System  Parameters.   General  flow  charts  describing  the 

logic  included  in  each  event  of  the  simulation  plus  the 
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interrelationship  of  events  are  included  as  Figures  9  through 
12.   Detailed  flow  charts  for  the  entire  model  are  provided 
in  Appendix  B.   Random  observation  is  the  subroutine  that 
generates  random  points  in  time  at  which  certain  parameter 
values  are  recorded  in  order  to  obtain  certain  useful  sta- 
tistics about  the  system.   Calculation  of  Measures  of 
Performance  is  the  subroutine  that,  at  the  completion  of 
any  unit  period  (a  time  period  which  represents  10  years 
of  actual  system  operations) ,  interrogates  the  information 
contained  in  all  system  parameters  and  transfers  control  of 
the  computer  program  to  the  calculation  phase  of  the  measures 
of  performance  listed  in  Table  I.   Incrementation  of  system 
parameters  is  the  subroutine  which  increments  the  current 
value  of  R,  the  reorder  point,  by  one.   It  also  clears  the 
results  of  the  last  unit  period  from  the  program  and  resets 
all  the  random  number  generators  to  the  same  values  used  in 
the  previous  run  to  ensure  meaningful  comparisons  between 
runs . 

C.   INPUT 

All  inputs  to  the  computer  program  are  contained  on 
three  punched  cards  called  Initial  Cards  and  four  Function 
Follower  Cards.   The  format  for  this  information  and  the 
units  to  be  used  are  detailed  in  Tables  XIV  through  XVI. 
With  one  set  of  input  cards  a  user  can  simulate  the  inter- 
actions involved  for  as  long  as  140  years  or  14  unit  periods 
of  ten  years  each.   Each  unit  period  will  increment  the 


54 


(8)  I    RECEIPT 

RECEIPT    D  OF  ITEMS 

.MATFER.ALi-R-E-P-AiR-E-? 


(6) 

BACKORDER 
B 


I 

S 

S 

u 


ON  HAND     (D-B) 


PROCUREMENT 
SOURCE 


Figure  9.    Relationships  of  events. 
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value  of  R  by  one  starting  with  the  initial  value  chosen  by 
the  user  (generally  zero) .   In  order  for  the  user  to  vary 
parameter  values  other  than  the  reorder  level,  a  new  set  of 
initial  cards  must  be  prepared  by  the  user  for  the  program. 
All  input  values  must  be  integers  and  in  units  of  hours  if 
it  is  a  time  input.   For  example  in  the  three  columns  allotted 
for  Q  on  card  one  the  user  could  enter  anything  from  001  to 
999  for  one  item  up  to  999  items  or  in  the  columns  allotted 
for^  ,  the  mean  time  between  demands  (failures)  ,  on  card 
two  the  user  could  enter  anything  from  999  to  001  to  repre- 
sent failures  rates  of  .7308  per  month  (one  every  999  hours) 
to  730  per  month  (one  every  hour) . 

D.   FORMAT  SECTION 

This  section  contains  the  tables  which  identify  and 
explain  the  various  system  parameters  and  characteristics 
needed  as  user  input  in  order  to  use  the  simulation.   It 
also  includes  the  feasible  range  for  each  entity.   This 
input  must  be  typed  by  the  user  on  three  Initial  Cards  and 
four  Function  Follower  Cards. 

1.  Description  Of  Entities  Found  On  The  Initial  Cards 
This  input  is  summarized  in  Tables  XIV  through  XVI. 

2 .  Description  Of  Function  Follower  Cards  One  and  Two 
These  cards  contain  the  function  values,  Y-,  and 

argument  values,  X.,  for  the  empirical  distribution  of 
repair  time.   This  is  a  discrete  numerical  valued  function 
wjLth  up  to  ten  different  points.   Figure  12  shows  that  this 
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TABLE  XIV.   FIRST  INITIAL  CARD 


Q       Reorder  quantity.   Number  of  items  reordered  at 
one  time. 

R       Reorder  Point.   If  Inventory  Position  is  <£,  R  an 
order  is  placed. 

L       Minimum  Batch  Size.   Under  the  batch  repair  policy, 
repair  won't  start  unless  at  least  L  units  are 
currently  waiting  to  be  repaired.   Ensure  L=0  if 
the  infinite  server  queueing  policy  for  repair  is 
desired. 

K       Maximum  Batch  Size.   Under  the  batch  repair  policy, 
this  is  the  maximum  number  of  items  that  can  be 
repaired  at  one  time. 

P       Probability  of  Repair.   The  probability  that  a 

failed  item  can  be  repaired.   If  P=0000,  (P=0.0) 
the  model  degenerates  to  the  consumables  case. 
If  P=1000,  (P=1.0)  the  model  becomes  one  without 
degredation  of  supply  which  implies  that  no  orders 
are  placed  and  no  outside  procurement  is  needed. 

.       Mean  time  between  failures.   The  expected  number  of 
^D  hours  between  demand. 

FLAG/R  Repair  Time  Distribution  Flag.   One  of  four  numbers, 

3  Constant  repair  time. 

4  Uniform  distribution  for  repair  time. 

5  Normal  distribution  for  repair  time. 

6  Empirical  distribution  for  repair  time. 
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TABLE  XV.   SECOND  INITIAL  CARD 


aj  Mean  Time  For  Repair.   The  expected  time  required 

R     for  repair  of  an  item  or  batch  of  items. 

LBR     Minimum  Repair  Time.   The  lower  bound  on  repair 
time  for  the  uniform  distribution  of  repair  time 
(used  when  FLAG/R=4) . 

UBR     Maximum  Repair  Time.   The  upper  bound  on  repair 

time  for  the  uniform  distribution  of  repair  time, 
(used  when  FLAG/R=4) . 

^  Standard  Deviation  For  Repair  Time.   The  standard 

Ufl    deviation  for  the  normal  distribution  of  repair 
time,  (used  when  FLAG/R=5) . 

FLAG/R  Lead  Time  Distribution  Flag.   One  of  four  numbers 

7  Constant  lead  time. 

8  Uniform  Distribution  for  lead  time. 

9  Normal  Distribution  for  lead  time. 

10  Empirical  Distribution  for  lead  time. 


^L 


Mean  Time  For  Procurement.   The  expected  value  for 
procurement  lead  time. 


TABLE  XVI.   THIRD  INITIAL  CARD 


LBL     Minimum  Lead  Time.   The  lower  bound  on  lead  time 
for  the  uniform  distribution  of  lead  time. 

UBL     Maximum  Lead  Time.   The  upper  bound  on  lead  time 
for  the  uniform  distribution  of  lead  time. 


6L 


Standard  Deviation  For  Lead  Time.   The  standard 
deviation  for  the  normal  distribution  of  lead  time. 


XH17    Type  Of  Output  Desired  Flag.   One  of  two  numbers. 

Zero  for  summary  output,  or  one  for  standard  output. 
(Normally  a  zero  should  be  used  to  conserve  printing 
time . ) 
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function  has  the  same  function  value,  FN.=Y  ,  for  all 

3      l 

argument  values  in  the  interval  between  X.-l  and  X..   No 

3        3 

interpolation  is  performed,  and  the  vlaue  at  the  right- 
hand  end  of  the  interval  is  used.   For  argument  values 
less  than  the  magnitude  of  X  ,  FN ,  =Y, .   For  argument  values 

greater  than  the  last  value  of  X.,  FN  =Y  . 
6  j  ■   n  n 

The  function  follower  cards  contain  the  individual 
points  (successive  pairs  of  X.,  Y.  values)  of  the  function. 
For  example,  consider  the  function  pictured  in  Figure  12  as 
the  repair  time  distribution,  the  function  follower  cards 
would  then  be:   0.06,439/0.07,439/0.08,439/0.58,464/0.66,486/ 
0.68,558/0.69,558/0.71,558/0.96,582/1.0,606. 

Notice  that  in  this  example  there  are  really  only- 
six  different  function  values  FN.,  so  argument  values 
.06,  .07,  .68,  and  .69  are  dummy  values  and  are  put  in  only 
to  meet  the  requirement  of  having  ten  pairs  of  X.,Y.  values. 
3.   Description  Of  Function  Follower  Cards  Three  and  Four 

These  cards  are  used  exactly  the  same  way  as  Function 
Follower  Cards  One  and  Two  except  they  contain  the  function 
and  argument  values  for  the  empirical  distribution  of  lead 
time.   For  example,  if  the  function  pictured  in  Figure  13 
were  the  lead  time  distribution,  FFC ' s  three  and  four  would 
be:   0.2,360/0.4,384/0.5,432/0.6,480/0.8,768/0.9,840/ 
0.95,864/0.99,888/0.995,912/0.999,984. 

Notice  that  once  FLAG/R  and  FLAG/L  have  been 
selected  the  program  ignores  those  cards  which  involve  the 
parameters  of  the  other  distributions.   Therefore  only  those 
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values  or  cards  which  affect  the  parameter  values  of  the 
particular  type  of  distribution  specified  need  be  worried 
about.   For  example,  if  FLAG/R=5  and  FLAG/L=10,  the  user 
must  ensure /J^,   01  andy^L  have  the  values  desired  and  FCC's 
three  and  four  have  the  values  desired.   The  values  of  LBR, 
UBR,  LBL,  UBL,  o"^  ,  and  FFC '  s  one  and  two  can  be  ignored. 
4 .   Examples  of  Initial  Cards 

Tables  XVII  through  XIX  give  examples  of  the  three 
Initial  Cards. 


TABLE  XVII.   EXAMPLE 

OF  INITIAL  CARD  ( 

3NE 

Columns 

Name 

Range 

Units 

8-14 

INITIAL 

19-22 
23-24 

XH1, 
Q 

1-99 

ITEMS 

25-29 
30-31 

/XH2, 
R 

0-99 

ITEMS 

32-36 
37-38 

/XH3, 
L 

0-K 

ITEMS 

39-43 
44-45 

/XH4, 

K 

0-99 

ITEMS 

46-50 
51-54 

/XH5, 
P 

0-1000 

.001 

55-59 
60-62 

/XH6, 

0-999 

HOURS 

63-67 
68-69 

/XH7, 
FLAG/R 

3,4,5 ,  or  6 

NONE 

All  input  information  is  right- justified  in  the 
space  allotted.   Example: 


COLUMN  NUMBER  - 


28 
t 


39 

t 


55 


69 
t 


INITIAL   XH1,10/XH2,01/XH3,00/XH4,01/XH5,0900/XH6,0  61/XH7,03 


64 


TABLE  XVIII.   EXAMPLE  OF  INITIAL  CARD  TWO 


Columns 

8-14 

19-22 
23-25 

26-30 
31-33 

34-39 
40-43 

44-49 
50-52 

53-58 
59-60 

61-66 
67-69 


Name 

INITIAL 

XH8, 

/XH9, 
LBR 

/XH10, 
UBR 

/XH11, 


/XH12, 
FLAG/L 

/XH13, 


Range 

0-999 

0-998 

0-998 

0-580 

7,8,9,  or  10 

0-999 


Units 

HOURS 

HOURS 

HOURS 

HOURS 

NONE 

HOURS 


All  input  information  is  right- justified  in  the 
space  allotted.   Example: 


COLUMN  NUMBER 


53 
t 
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INITIAL   XH8,336/SH9,lt38/XH10,0504/XHll,095/XH12,07/XH13,504 


TABLE  XIX.   EXAMPLE  OF  INITIAL  CARD  THREE 


Columns 

Name 

8-14 

INITIAL 

19-23 
24-26 

XH14, 
LBL 

27-32 
33-36 

/XH15, 
UBL 

37-42 
43-45 

/XH16, 

46-51 
52 

/XH17, 

Range 

0-998 
0-1998 
0-580 
0,  or  1 


Units 

HOURS 
HOURS 
HOURS 
NONE 


All  input  information  is  right- justified  in  the 
space  allotted.   Example: 

COLUMN  NUMBER  -  26  37  51 

t  *  t 

INITIAL    XH14,240/XH15,0  768/XH16,116/XH17,0 
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E .   OUTPUT 

The  computer  program  has  two  different  sets  of  output, 
summary  (special)  and  standard.   The  type  of  output  received 
depends  on  the  value  of  XH17.   If  XH17=0  the  summary  output 
is  received,  but  if  XH17=1  the  standard  output  is  received. 
Computer  running  and  printing  time  is  greatly  increased  if 
the  standard  output  is  used.   Therefore,  it  is  recommended 
that  XH17  be  set  to  zero.   If  there  is  a  major  error  in  the 
program,  the  standard  output  will  be  printed  out  automatic- 
ally.  Following  is  an  explanation  of  the  two  sets  of  output 

1.   Summary  Output 

An  example  of  the  summary  output  is  illustrated  in 
Appendix  D.   Notice  that  there  are  three  major  divisions: 
(1)  savevalues,  (2)  storage  statistics  and  (3)  queue 
statistics . 

a.   Savevalues 

There  are  four  subdivisions  in  this  part.   The 
first  is  a  list  of  nonzero  halfword  savevalues  between 
five  and  100  for  system  parameters  and  statistics.   This 
list  is  included  for  completeness  only  and  not  because  of 
their  nature  or  importance.   The  second  subdivision,  Clock 
Time,  merely  shows  the  number  of  hours  of  operation 
simulated.   Normally  this  will  be  87600  hours  or  ten  years 
of  operation,  which  has  been  defined  to  be  a  unit  period. 
The  third  subdivision,  Critical  System  Parameters,  lists 
the  values  of  Q,  K,  and  L  which  the  user  has  chosen  and 
input  on  the  initial  cards.   It  also  prints  out  the  value 
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of  R  for  this  particular  time  period.   The  last  subdivision, 
Measures  of  Performance,  includes  most  of  the  important 
statistics  to  be  considered  when  judging  or  comparing 
alternative  inventory  systems. 
b.   Storage  Statistics 

This  division  provides  the  important  statistics 
for  the  order  department  (ORDER) ,  the  infinite  server 
queueing  repair  process  (FASTR) ,  and  the  single  server 
queueing  repair  process  (SLOWR) .   The  important  items  and 
explanation  of  each  storage  follow.   (1)  Storage  --  the  name 
used  to  represent  the  physical  entity.   (2)  Capacity  --  the 
maximum  number  of  transactions  which  can  use  the  storage 
at  any  one  time.   (3)  Average  Contents  --  the  average 
number  of  transactions  that  were  in  the  storage  during  the 
run.   (4)  Entries  --  the  number  of  transactions  which 
entered  the  storage.   For  SLOWR  this  is  the  number  of 
batches  of  items  that  entered  the  repair  process  rather 
than  the  number  that  were  repaired  as  in  the  FASTR  storage. 
Likewise,  it  is  the  number  of  orders  placed  rather  than  the 
number  of  items  ordered  for  the  ORDER  storage.   (5)  Average 
Time/Tran  --  the  average  number  of  hours  that  an  order  was 
outstanding,  or  the  average  number  of  hours  that  it  took 
for  an  item  to  be  repaired,  or  the  average  number  of  hours 
it  took  for  a  batch  of  items  to  be  repaired.   (6)  Maximum 
Contents  --  the  maximum  number  of  transactions  that  used 
the  storage  at  any  time  during  this  run. 
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c.   Queue  Statistics 

This  division  reflects  the  important  statistics 
for  the  infinite  server  queue  repair  department  (FSTRQ)  and 
the  single  server  queue  repair  department  (SLOWQ) .   The 
important  items  and  explanation  of  each  follow.   (1)  Queue  -■ 
the  name  used  to  represent  the  physical  entity.   (2)  The 
explanation  for  maximum  contents,  average  contents,  and 
average  time/tran  are  the  same  as  for  storages.   (3)  Total 
Entries  --  the  number  of  transactions  which  entered  the 
queue.   (4)  Zero  Entries  --  the  number  of  transactions 
which  entered  the  queue  and  departed  with  no  time  delay. 
Note  that  for  the  FSTRQ  it  should  always  be  the  same  as  the 
total  entries.   (5)  Average  Time/Trans  --  the  average  number 
of  hours  that  the  delayed  transactions,  only,  spent  waiting 
to  be  repaired. 

2 .   Standard  Output 

An  example  of  the  standard  output  is  illustrated  in 
Appendix  D.   Notice  that  the  standard  output  includes  all  of 
the  items  listed  in  the  summary  output  as  well  as  the  follow- 
ing items.   (1)  table  statistics,  (2)  group  counts,  and 
(3)  block  counts . 

The  table  statistics  provide  the  statistics  of  the 
following  tables:   BTCHZ ,  BTIME,  ENOH,  NBARB,  NUMC ,  and 
RDYRT.   The  important  items  and  explanations  of  each  follow. 
(1)  Table  --  the  name  used  to  collect  the  particular 
statistic.   (2)  Entries  in  Table  --  the  number  of  times 
that  an  entry  was  made  in  the  table.   (3)  Mean  Argument  -- 
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the  average  value  of  the  entries  made  in  the  table. 
(4)  Standard  Deviation  --  the  standard  deviation  of  the 
entries  made  in  the  table.   (5)  Upper  Limit  --  the  maximum 
value  of  each  cell  of  the  table.   For  example,  all  entries 
between  51  and  100  fall  in  cell  two  of  BTIME  table. 
(6)  Observed  Frequency  --  the  number  of  entries  which  fall 
into  this  cell.   (7)  Percentage  of  Total  and  Cumulative 
Percentage  are  self-explanatory.   (8)  Multiple  of  Mean  -- 
the  upper  limit  of  the  cell  expressed  as  a  multiple  of  the 
mean  of  the  table.   The  additional  pieces  of  information 
given  in  the  tables  as  well  as  the  queue  and  storage 
statistics  which  are  not  explained  in  detail  here  are  those 
which,  under  normal  circumstances,  are  of  no  concern  to  the 
user.   However,  there  may  be  occasions  when  a  user  wishes 
to  make  changes  internal  to  the  program  such  as  adding  the 
capability  of  handling  compound  Poisson  demands.   In  such 
a  case  the  standard  output  would  be  a  necessity,  especially 
the  block  counts,  in  helping  the  user  debug  the  program. 
The  reason  this  output  was  included  is  to  show  the  user 
that  this  additional  information  is  available  from  the 
computer,  and  to  show  the  user  the  type  of  output  he  will 
receive  by  default  should  an  error  be  made  while  using  the 
model.   If  the  user  desires  to  receive  the  standard  output, 
XH17  should  be  set  to  one  on  initial  card  three.   If  one 
desires  further  information  relating  to  computer  output, 
he  should  see  the  IBM  User's  Manual  [Ref.  4]. 
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VI.   USES  AND  APPLICATIONS  OF  THE  MODEL 

This  model  was  built  for  three  major  reasons.   First 
and  foremost  was  to  give  a  convenient  and  useful  tool  to 
present  day  inventory  managers  faced  with  the  decisions  of 
how  much  to  buy,  when  to  buy,  and  how  to  conduct  repairs 
for  the  items  under  their  cognizance.   As  was  pointed  out 
earlier,  there  are  presently  two  major  approaches  in  vogue 
to  help  the  manager  solve  his  problem;  cost  analysis  and  the 
use  of  simple  procurement  policies  with  specified  measures 
of  supply  performance.   Under  the  cost  analysis  approach 
the  user  runs  into  the  difficult  problems  of  quantifying 
costs.   For  instance  I,  the  inventory  holding  cost,  consists 
of  many  things,  some  of  which  all  managers  agree  on,  such 
as  storage  and  handling,  and  some  which  they  don't  agree 
on,  such  as  loss,  breakage,  and  opportunity  cost.   Similar 
disagreements  concern  the  fixed  order  cost  as  well.   For 
example,  the  Navy  uses  a  sliding  scale  from  $25  to  $75, 
while  the  Army  uses  a  different  sliding  scale  from  $50  to 
$400.   Further  problems  arise  when  one  tries  to  assign  a 
cost  to  a  critical  item  being  out  of  stock  when  needed. 
Primarily  because  of  these  difficulties  in  quantifying 
costs,  the  second  approach  was  used  in  this  paper.   It 
was  felt  that  statements  such  as  "a  fill  rate  of  95%  must 
be  maintained"  would  at  least  mean  the  same  thing  to  all 
managers  as  well  as  being  something  that  all  managers  could 
agree  on. 
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The  second  reason  this  model  was  built  was  to  provide 
an  efficient  way  of  obtaining  insight  into  the  problems 
faced  by  the  inventory  manager.   Simulation  is  a  straight- 
forward and  logical  way  of  doing  this.   In  addition  simula- 
tion avoids  some  of  the  pitfalls  of  the  more  classical 
means  of  analysis.   For  example,  the  model  can  use  any 
time-step  desired  and  simulate  as  far  into  the  future  as 
desired.   The  model  easily  handled  the  different  types  of 
repair  processes  which  analytically  prove  to  be  quite 
cumbersome.   When  these  reasons  are  coupled  with  the 
combined  effects  of  uncertainty,  dynamic  interaction  of 
system  events,  and  interdependence  among  system  parameters, 
it  is  easily  seen  why  this  problem  was  particularly  amenable 
to  simulation.   Richards  [Ref.  7]  showed  that  analytical 
expressions  could  be  obtained  for  the  individual  measures 
of  performance  used  in  this  paper.   However,  to  actually 
generate  numbers  from  these  expressions  for  a  sensitivity 
analysis  of  this  size  would  be  a  monumental  task.   Therefore, 
simulation  again  seemed  to  be  called  for.   Also,  the  advan- 
tage of  time  compression,  both  in  being  able  to  simulate 
into  the  future  in  a  matter  of  seconds  and  being  able  to 
obtain  numerical  values  for  the  various  measures  of  per- 
formance, made  simulation  appealing. 

A  major  criticism  of  simulation  is  that  the  simulated 
queueing  model  yields  only  an  estimate  of  parameters, 
such  as  a  waiting  line's  average  length  or  associated 
probability  of  delay,  and  as  such  is  subject  to  statistical 
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error.   A  second  criticism  is  that  if  the  system  is  so 
complicated  and  complex  that  the  other  methods  fail  and  a 
simulation  is  required,  it  will  necessarily  be  a  very 
difficult  job  to  build  the  model  and  analyze  its  results. 

In  this  particular  case  the  first  disadvantage  acted 
to  a  much  lesser  degree  since  analytical  expressions  for 
all  the  measures  of  performance  used  in  this  model  do  exist. 
This  fact  allowed  less  replication  to  obtain  trustworthy 
answers  in  the  search  for  the  best  choices  of  the  policy 
parameters.   A  moderate  number  of  runs  in  a  reasonable 
period  of  time  accompanied  with  the  statistical  tests  of 
section  four  also  lent  credibility  to  the  answers  and 
conclusions  obtained.   The  second  disadvantage  noted  went 
for  naught  since  the  model  exists  and  works. 

Finally,  the  third  reason  that  this  model  was  constructed 
was  to  demonstrate  the  validity  of  certain  analytical 
expressions  for  some  of  the  measures  of  performance  used 
in  this  paper.   As  was  pointed  out  in  section  four,  the 
model  succeeded  in  doing  this.   From  the  comments  above 
and  the  results  presented  in  this  paper,  it  is  easily  seen 
that  the  model  can  be  used  effectively  to  describe  current 
inventory  systems,  to  explore  hypothetical  inventory 
systems,  and  to  design  or  improve  given  inventory  systems, 
as  well  as  be  used  as  a  teaching  aid.   The  simulation 
model  can  be  used  to  help  answer  such  a  current  question 
as  how  to  best  conduct  repairs  on  the  engines  for  the  new 
class  of  destroyers  DD963. 
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In  applying  this  model  to  any  of  the  uses  mentioned 
above,  one  is  always  faced  with  the  problem  of  deciding 
which  measure  of  performance  to  use  as  the  criterion  for 
making  decisions.   In  this  paper  the  expected  number  of 
unit  years  of  backorders  per  unit  period  was  used.   One 
reason  it  was  used  over  the  average  number  of  items  on 
backorder  per  unit  period  was  that  it  considers  both  the 
number  of  backorders  and  the  amount  of  time  the  backorders 
were  outstanding,  while  the  other  merely  considers  the 
number  of  backorders,  totally  ignoring  the  amount  of  time 
the  backorders  were  outstanding.   It  was  felt  that  both 
the  fill  rate  and  the  ready  rate  were  even  worse  than  the 
average  number  of  backorders  since  they  don't  even  take 
into  consideration  the  number  of  units  backordered,  but 
merely  note  whether  or  not  any  backorders  exist.   Perhaps 
the  most  significant  reason  the  expected  number  of  unit 
years  of  backorders  per  unit  period  was  used  as  the  criterion 
was  that  once  it  was  considered  acceptable,  the  remaining 
measures  of  performance  were  also  acceptable. 

There  were  occasions  when  another  measure  of  performance, 
the  average  number  of  items  on  hand  at  an  arbitrary  point 
in  time,  was  used  to  help  determine  an  optimal  operating 
point  for  a  given  system.   Such  was  the  case  in  experiment 
one  when  there  appeared  to  be  a  number  of  acceptable 
operating  points,  (Q,R)  pairs,  for  the  system  under  investi- 
gation.  This  measure  of  performance  was  used  in  these 
cTises  to  differentiate  between  these  pairs. 
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As  mentioned  earlier,  this  model  can  also  be  used  as  a 
teaching  aid  as  well  as  demonstrating  analytical  results 
obtained  by  others.   Table  XX  (page  75)  is  a  good  example 
of  this.   If  one  sets  the  probability  of  repair  equal  to 
one,  the  model  degenerates  to  the  all  repairables  case 
without  degradation.   From  the  table,  one  can  see  how  the 
various  measures  of  performance  change  as  the  mean  demand 
during  repair  time  varies.   Further,  if  Q  is  required  to 
be  one,  the  results  of  this  model  can  be  compared  with  the 
results  of  Feeney  and  Sherbrooke  [Ref.  6].   Notice  that  in 
their  results  the  fill  rate  was  always  zero  when  the  reorder 
point  was  zero.   This  result  was  a  consequence  of  their 
initial  conditions,  whereas  the  results  obtained  by  this 
model  do  not  include  that  restriction.   If  we  again  restrict 
the  model  by  considering  the  all  consumables  case  (i.e. 
letting  the  probability  of  repair  be  equal  to  zero)  the 
model  becomes  that  described  by  Hadley  and  Whitten [Ref .  3], 
Using  a  sample  problem,  one  again  has  an  example  of  how  the 
model  can  be  used  to  verify  known  analytical  results.   This 
example  and  comparison  to  the  cost  approach  is  summarized 

in  Table  XXI. 

Of  course  there  are  endless  experiments  one  could  carry 
out  with  the  use  of  this  model.   For  instance  one  could  use 
this  model  to  help  decide  which  items  should  be  designated 
repairable.   This  could  be  done  by  comparing  the  measures 
of  performance  obtained  by  considering  an  item  first  as  a 
consumable  and  then  as  a  repairable.   It  would  also  be  easy 
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to  extend  this  model  to  make  use  of  the  variance  reduction 
techniques  explained  in  reference  nine  or  to  add  the 
capability  of  handling  compound  Poisson  demands. 

VII.   SUMMARY  AND  CONCLUSION 

The  analysis  of  the  repairable  inventory  problem  pre- 
sented in  this  paper  was  based  upon  a  simulation  model 
thought  to  be  sufficiently  realistic  and  general  enough  to 
represent  many  different  inventory  systems. 

The  results  have  shown  that  the  simulation  can  be  a 
helpful  tool  for  the  managers  of  inventory  systems  to  use 
to  make  wise  decisions  governing  their  policies  concerning 
when  and  how  much  to  reorder.   It  was  demonstrated  in  this 
paper  that  the  mean  time  to  repair  an  item,  but  not  the 
repair  time  distribution  itself,  was  extremely  critical  to 
system  performance.   It  was  demonstrated,  at  least  for 
systems  with  a  mean  repair  time  less  than  or  equal  to  1.4, 
that  a  constraint  placed  on  the  minimum  repair  batch  size 
had  a  negligible  effect  on  system  performance.   It  was 
further  demonstrated  for  this  same  system  that  reducing 
the  mean  repair  time  by  40%  was  more  profitable  than  doubling 
the  maximum  repair  batch  size.   It  also  demonstrated  that 
increasing  maximum  repair  batch  size  beyond  those  values 
which  gave  a  service  rate  of  .417  or  less  had  negligible 
effect  in  terms  of  increasing  system  performance.   This 
paper  also  emphasized  the  power  of  using  the  expected 
number  of  unit  years  of  backorders  per  unit  period  as  the 
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criterion  for  judging  system  performance.   This  measure, 
when  reduced  to  an  acceptable  level,  was  accompanied  by 
good  values  of  all  other  measures  of  performance  as  well. 
This  is  not  true  of  the  other  measures.   For  example,  if 
the  fill  rate  had  been  used,  Table  XIII  shows  that  an 
acceptable  fill  rate  doesn't  necessarily  imply  that  all 
the  other  measures  of  performance  are  also  acceptable. 

It  is  believed  that  this  paper  supplements  the  tools 
which  are  needed  to  make  intelligent  decisions  concerning 
the  repairable  item  inventory  system.   Of  course  it  can 
also  be  used  for  the  consumables  case  by  merely  setting  the 
probability  of  repair  equal  to  zero. 

A  couple  of  further  remarks  seem  to  be  in  order.   First, 
it  would  seem  worthwhile  to  determine  if  current  repair 
departments  could  be  utilized  more  fully.   Perhaps  this 
could  be  achieved  by  spending  a  greater  fraction  of  budget 
dollars  to  repair  available  carcasses  rather  than  continuing 
to  procure  new  material.   Secondly,  this  study  seems  to 
indicate  that  study  of  a  stock  point's  physical  plant  for 
repairables  could  be  profitable.   For  example,  the  fourth 
experiment  indicated  that  the  minimum  repair  batch  size 
required  before  repair  begins  on  a  batch  of  items  could 
fluctuate  to  some  degree  before  system  performance  is 
degredated. 

Much  has  been  learned  about  the  interrelationships  and 
effects  of  many  of  the  parameters  of  repairable  and  non- 
repairable  inventory  systems.   Perhaps  this  has  been  the 
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most  significant  result  of  the  project.   It  would  seem 
impossible  for  anyone  who  would  work  with  this  simulation 
not  to  come  away  with  a  better  understanding  of  the  hidden 
mechanisms  and  inner  workings  of  an  inventory  system. 

In  closing,  a  simulation  model  for  continuous  review 
inventory  systems  was  presented.  The  model  used,  demonstrat 
ed  the  sensitivity  of  Q,  R,  K,  and  L  to  numerous  operational 
constraints  and  system  parameters.   A  class  of  simple 
policies  with  various  measures  of  performance  was  used 
rather  than  a  traditional  cost  postulated  model  to  determine 
acceptable  parameters  for  a  given  inventory  system.   It  is 
felt  that  the  model  can  be  successfully  implemented  in  a 
very  large  number  of  inventory  problems. 
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APPENDIX  A.   STATISTICAL  TESTS 

Analysis  of  variance.  This  test  was  used  to  find  the 
significant  difference  in  observed  values  of  the  expected 
number  of  unit  years  of  backorders  per  unit  period. 


Sum  of 
Squares 

Degrees 
Freedom 

of 

Mean 
Square 

F  Ratio 

Between  groups 

436.3584 

11 

39.6689 

194.0706 

Within  groups 

18.8052 

92 

0.2044 

Total 

455.1636 

103 

Table  value  from'  F-Distribution  for  92  and  11  degrees  of 
freedom  at  .01  significance  level  =  2.48. 

194.0706  >>   2.48,  therefore  the  null  hypothesis  was 
rejected. 

Null  hypothesis:  The  observed  values  for  the  expected 
number  of  unit  years  of  backorders  per  unit  period  are  the 
same  for  each  treatment  (reorder  point  value) . 
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K.  A.  Scheffe's  Multiple  Comparison  Test.  This  test 
was  used  to  find  the  magnitude  required  for  the  expected 
number  of  unit  years  of  backorders  per  unit  period  to  be 
significantly  different  from  zero. 

If      c  Ic.l  )2 

1  J' >.      F 

~  ^>     (l-o<)(n  )(n2) 

(N-l)O"  (C.)  )  L 

then  the  hypothesis  that  the  contrast  C.  differs  from 

zero  is  rejected. 

t 
c.   =      y         C.  .  T.  . 
J     i=l     1J   1 

A 

W  (C.)    is  the  estimated  variance  of  contrast  C. 

w  (T.)    is  the  estimated  variance  of  the  i   treatment=S  . 

T.        is  the  sum  of  the  n.  observations  for  the  i 
1  1 

treatment. 

2 
S         is  the  mean  square  of  the  experimental  error. 

t         is  the  number  of  treatments  included  in  contrast  C. 

C  .       are  the  coefficients  for  the  j    comparison  with 


ij 


the  i™  treatment 


c<        is  the  significance  level  chosen=.01  for  this 

particular  test. 
N-l       is  the  total  number  of  treatments  less  one. 
(n, ) (n  )   are  the  degrees  of  freedom. 
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Let  C  be  the  contrast  between  (QR)  pairs  (10,11)  and 
(10,10)  respectively.   Then: 


C1=C10  1T10+CH  1Tir^3H-25.3176)  +  (-4)(-22.1442)  =  12.62 

V(CJ=5"   C   2  (T.)  =  (32)  (8)  (.2044)  +  (-4)2(6)  (.  2044)  =34  .  339 
1   i=io   il     1 

flCiJ  )2       =     (12. 62)2      =   .422 
(N-l)tC(ci))      (11)(34.339) 

.422  J>       2.48  therefore,  H  is  not  rejected. 

Let  C-  be  the  contrast  between  (QR)  pairs  (10,11)  and 
(10,9)  respectively.   Then: 

C2=Cg  2T9+C11  2T11=(4)(22.404)+(-6)(-22.1442)=43.24 

V(C2)  =  (42)  (9)  (.2044)  + (-6)  2  (6)  (.  2044)  =  73  .  584 

(    IcJ    )2  =  (43. 24)2  =      2.309        2.309  j>   2.48 

(N-l) (<f(C2))  (11)  (73.584)  Therefore,    HQ 

is  not 
rejected. 

Let  C3  be  the  contrast  between  (QR)  pairs  (10,11)  and 

(10,8)  respectively.   Then: 

C3=C8  3T8+C11  3T11=(4)  ("14-353)  +  (-6)  (-22 . 1442) =75  .  453 
V  (C3)  =  (4)2(9) (.2044) +  (-6) 2  (6) (. 2044) =73 . 584 

CK|)2  2 

~  =  (75.455)z  =    6.994  6.994>   2.48 

(N-l) (^  (C3))  (11)  (73.584)  Therefore, 

reject  HQ ,  the 
null  hypothesis 
that  the  means 
are  equal. 
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Let  C  be  the  contrast  between  (QR)  pairs  (10,11)  and 
(10, m)  respectively.   Then: 

(  lCml  )2  ,   ,  2 

'  =  2.48  =    (  Cm  )Z  =  (2.48)(11)(73.584)  = 

(N-l)(A(Cm)) 

2007.372 


C   =  44.8   (minimum  contrast  value 

required  for  statistical 
significance.) 

Then:   44.  *-Cx%mTx+Cn  ^T^-  (4)TX+  (-6)  (22  .  404) 

Tx  =  22.016 

loge  BT  =  -2.450 

B  =  .086  the  significant  difference  of  the  expected  number 
of  unit  years  of  backorders  from  .01 

BT  =  . 01+. 086=0. 096  the  minimum  significant  difference 
from  zero 
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APPENDIX  B    Detailed  Flow  Charts 
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